Two-dimensional solitons on the surface of magnetic fluids 
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We report an observation of a stable soliton-like structure on the surface of a ferrofluid, gener- 
ated by a local perturbation in the hysteretic regime of the Rosensweig instability. Unlike other 
pattern-forming systems with localized 2D structures, magnetic fluids are characterized by energy 
conservation; hence their mechanism of soliton stabilization is different from the previously discussed 
gain/loss balance mechanism. The radioscopic measurements of the soliton's surface profile suggest 
that locking on the underlying periodic structure is instrumental in its stabilization. 

PACS numbers: 47.65.+a, 47.62.+q ,47.60.+i 



To date, stable solitary waves have been experimen- 
tally observed in a variety of one-dimensional and quasi- 
one-dimensional physical systems. In 2D, dispersive non- 
linear systems are prone to collapse instabilities and 
hence the 2D solitons turned out to be more elusive. 
(Here we use the term "soliton" in a broad physical sense, 
as a synonym of localized structure.) So far, the list of 
experimentally detectable 2D localized objects was con- 
fined mostly to vortices in superfluids, superconductors, 
and other media on one hand, and dissipative solitons in 
nonequilibrium systems on the other. While the stability 
of the former is due to their nontrivial topology, the latter 
come into being via the balance of strong dissipation and 
energy gain. Examples include current filaments in gas 
discharge systems 1]; oscillons in fluids and granular ma- 
terials 0; breathing spots in chemical reactions Q and 
feedback and cavity solitons in optics 4]. Despite some 
encouraging theoretical insights, the question of whether 
2D non topological solitons can arise in conservative sys- 
tems has remained open. 

In this Letter we report an experimental observation 
of a strongly localized, stable stationary soliton on the 
surface of magnetic fluid (MF) in a stationary magnetic 
field. MF is a dispersion of magnetic nanoparticles, and 
thus has a high relative permeability /i r p. This is a 
lossless system; a horizontal layer of MF in a vertically 
applied magnetic induction B is characterized by the en- 
ergy density 0, Q : 

F[h(x,y)} = ^-h 2 (x,y) - J dzB ^ r - 1 H MF (x,y,z) 

+o-^/l + (d x h(x,v)) 2 + (dyHx,y)) 2 - (1) 

Here p and a are the density and surface tension of 
the MF, h(x,y) the local height of the liquid layer, and 
^mf(^,1/,^) is the magnetic field in the presence of the 
MF. The three terms in Eq. represent the hydrostatic, 
magnetic and surface energy, respectively. As the surface 
profile deviates from the flat reference state, the first and 
last term grow whereas the magnetic energy decreases. 
For sufficiently large 5, thisjgives rise to the normal field, 
or Rosensweig, instability d, 

Our experimental setup is sketched in Fig.[I] A 



Teflon® vessel with the radius R = 60 mm and depth of 3 
mm |9| is filled with MF up to the brim and placed on the 
common axis midway between two Helmholtz coils. An 
x-ray tube is mounted above the center of the vessel at a 
distance of 1606 mm. The radiation transmitted through 
the fluid layer and the bottom of the vessel is recorded by 
an x-ray sensitive photodiode array detector (16 bit) con- 
nected to a computer. The full surface relief is then re- 
constructed from the calibrated radioscopic images. For 
details see [lOj. The experiments were performed with 
the magnetic fluid EMG901, Lot F121901 AX from Fer- 
rotec. Its material parameters have been measured to 
/i r 3.2, p 1.406gcm- 3 , and a 25±0.7mN/m. 
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FIG. 1: Sketch of the experimental setup. See text for details. 

Starting from a flat layer at B = OmT, we increase the 
external induction adiabatically in steps of 15 //T, paus- 
ing for 30s after each increase. As shown in Fig.El(a), 
a deformation of the surface of the liquid appears first 
at the edge of the vessel. This is due to the magnetic 
field gradient induced by the discontinuous magnetiza- 
tion at the edge of the liquid layer. Increasing the induc- 
tion further gives rise to a fully developed pattern of the 
Rosensweig instability, as shown in (b,c). 

We measured the top-to-bottom height A of the sta- 
tionary fluid pattern arising in the adiabatic increase 
and decrease of B. To avoid the edge-induced imperfec- 
tions in the character of the bifurcation, we only consider 
spikes located within 11 mm from the center of the dish. 
FigurelHl displays results obtained for 400 values of B. As 




FIG. 2: Surface reliefs as reconstructed from the radioscopic 
images for a) B = 8.922 mT, and b), c) 10.407 mT. Each color 
indicates a layer thickness of 1 mm. 



B is increased, a sudden transition to the upper branch 
occurs at B c = 9.025mT. For B > B c , the entire surface 
is covered by a lattice of liquid spikes, which is hexag- 
onal away from the boundary. Decreasing £?, the order 
parameter A remains on the upper branch all the way 
to B* = 8.076mT where it drops to the fiat reference 
level. Thus the diagram shows a subcritical bifurcation 
to hexagons. The solid and dashed lines display a fit 
to the roots of the corresponding amplitude equation of 
Ref. 0. 

To study the stability of the flat surface to local per- 
turbations (in the hysteretic regime), a small air coil with 
the inner diameter of 8 mm was placed under the center 
of the vessel (see Fig.[Q. This allows to increase, locally, 
the magnetic induction. A local pulse of B+ = 0.68mT 
added to the uniform field of B = 8.91mT, produces a 
single stationary spike of fluid, surrounded by a circular 
dip, which does not disperse after B+ has been turned 
off. Figure^] presents a measured relief of this radially- 
symmetric state which will be referred to as the soliton. 
The soliton is a stable nondecaying structure; it remained 
intact for days. After its formation at the center of the 
dish, the soliton was often seen to float around (with 
v ~ O.lmm/s), until reaching an equilibrium position 
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FIG. 3: The amplitude of the pattern for r < 11mm versus 
the magnetic induction. The crosses (dots) mark the val- 
ues for increasing (decreasing) induction respectively. The 
solid (dashed) lines display the least square fit to the roots 
A± = [7(1 + e) ± V7 2 (l + e) 2 + 4eg]/(2g), of the amplitude 
equation eA + 7(1 + e)A 2 - gA 3 = of Ref . Q with 7 = 0.281 
and g = 0.062. The full circles (squares) give the amplitude 
of the localized spike initiated at B = 8.91mT for increasing 
(decreasing) induction respectively. 




FIG. 4: A single soliton surrounded by the unperturbed 
magnetic liquid. The magnetic induction generated by the 
Helmholtz coils amounts to B = 8.91mT. The amplitude of 
the local pulse which produced the soliton was B+ = 0.68mT 
at the bottom of the vessel. 



somewhere near the edge of the dish. This behavior can 
be attributed to radial gradients of the magnetic field due 
to the discontinuous magnetization at the edge and the 
ring of spikes pinned along the perimeter of the dish. 

We examined the range of stability of the soliton gen- 
erated by a pulse with B+ = 0.68 mT added to the 
uniform induction B — 8.91mT. Reducing B adiabat- 
ically we measured the corresponding amplitude of the 
soliton (marked by full squares in Fig.[3J). Similarly to 
the spikes in the hexagonal pattern, the height of the 
soliton decreases as B is reduced. The soliton decays 
for B < 8.09mT, which is close to B* = 8.076 mT, the 
lower stability boundary of the hexagonal pattern. As 
B is increased, the amplitude of the soliton grows, as 
indicated in Fig.[3| by full circles. At B = 9.055mT, 
a sudden transition from the soliton to the fully devel- 
oped Rosensweig pattern occurs. This value is somewhat 
larger than B c = 9.025 mT; this is due to the fact that 
the birth of the soliton produces a slight drop in the fiat 
layer thickness, which shifts B c - similarly to Fig. 5 in 
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FIG. 5: The filled squares mark the profile of one period of 
the hexagonal pattern, measured at B — 9.07mT in the center 
of the vessel; r < 8.8mm. Azimuthally averaged height pro- 
files of two different solitons, measured at the same induction 
are depicted by open symbols (one) and a dashed line (the 
other) . 

Ref.0. 

In order to illustrate the robustness of the soliton's 
shape we show in Fig.EI the azimuthally averaged pro- 
files of two different solitons, produced in two separate 
experiments at B = 9.07mT . The profiles are practically 
indistinguishable. Also plotted are two half-periods of 
the corresponding hexagonal lattice. In agreement with 
Fig.0 the soliton is about 1mm taller than the spikes 
of the lattice. This may be attributed to the fact that 
the spikes emerge simultaneously, and thus have to share 
the liquid available. However, the width of the soliton 
is exactly equal to the period of the lattice. Therefore, 
there is a preferred wavelength in the system, defined by 
the lattice, to which the soliton locks. As we show now, 
this width locking is the central part of the soliton's sta- 
bilization mechanism. 

Consider the dispersion relation of a semi-infinite layer 
of inviscid MF 0, 

u? = gk - /i /i r (/ir ~ l)2 -H 2 k 2 + -k\ (2) 

fJL r + 1 p p 

Here uo is the frequency, k = \k\ the wavenumber, H 
the strength of the external magnetic field, and po the 
magnetic field constant. While the first and the third 
terms account for the gravity and capillary effects and 
are common for all fluids, the second term is specific just 
for the MF. As H is increased above H C1 where H 2 = 
(p9 cr ) 1 ^ 2 \ a band of wavenumbers with uj 2 < 

appears around k c = (pg/a) 1 ^ 2 and the flat state loses 
its stability to the hexagonal pattern. Since the soliton 
should decay to the fiat surface as x 2 — > oo, there can be 
no stable solitons for H > H c . 

Let now H < H c . On the qualitative level, our system 
can be modelled by a conservative analogue of the Swift- 
Hohenberg equation: 

u + (fcg + V 2 ) 2 ^ + au = 3bu 2 - 2cu 3 ; 6, c> 0. (3) 



Eq.((3J) can be used as a model since (a) it has a nonmono- 
tonic dispersion relation uo 2 = a + k$ — 2k^k 2 + /c 4 , where, 
as in (2), the destabilizing /c 2 -term is opposed by a higher 
power of k; and (b) it has a symmetry-breaking hys- 
teretic nonlinearity which was shown to provide a fairly 
accurate approximation of the amplitude of the hexag- 
onal pattern 7] (see FigO) We have verified, numer- 
ically, that Eq.© does indeed have a stable stationary 
radially-symmetric soliton solution coexisting with stable 
hexagons in a broad parameter range. Its stability can 
be explained by a Derrick-type argument for the corre- 
sponding energy functional, 

E = J[ii 2 + V(u) - 2k 2 (Vu) 2 + (V 2 u) 2 ]d 2 x. (4) 

Here V(u) = (a + kfyu 2 — 2bu 3 + cu 4 . A scaling pertur- 
bation u(x) — > u(kx) takes the stationary energy to 

E(n 2 ) = ^2 J Vd 2 x-2kl J (Vu) 2 d 2 x+n 2 J (V 2 u) 2 d 2 x. 

The first term (nonlinearity) opposes the dispersive 
broadening of the soliton (for which k — > 0) while the last 
one prevents the nonlinear blow-up (for which k — > oo). 
In a similar way, the first and last term in J2J make 
contributions to the energy which scale as n~ x and ft, 
respectively. The first term (along with the nonlinear- 
ity) opposes the spreading and the last one arrests the 
blow-up. Next, setting the derivative (dE / 'dn 2 ) K= i to 
zero, gives J Vd 2 x = J (V 2 u) 2 d 2 x. Using this relation, 
the second derivative, d 2 E/d(K 2 ) 2 , is calculated to be 
2 j(y 2 u) 2 d 2 x > 0, which means that the soliton ren- 
ders the energy minimum. If Eq.© did not include the 
higher-derivative term, the energy would not have a non- 
trivial minimum. The introduction of a higher derivative 
(or, equivalently, the nonmonotonicity of the dispersion 
curve) sets a preferential wavelength in the system — to 
which the soliton locks and stabilizes. A similar stabi- 
lization mechanism was discussed before in the context 
of the wave front locking 11]; see also [T^ . 




FIG. 6: Nine solitons at B = 8.91mT. 



Applying, repeatedly, pulses of B+ and allowing the 
newly born solitons to drift away from the site of the 
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probe coil, we were able to generate two, three, and more 
solitons. Figure El presents an example of a 9-soliton con- 
figuration, with only one remaining at the center. In this 
way, it is possible to increase the surface energy of the 
liquid layer in steps. This is illustrated in Fig.0 which 
also shows the surface energy of the Rosensweig pattern 
as a hysteretic function of B. Thus, one can reach the 
region between the two branches of this function which 
is not accessible for the standard Rosensweig instability. 




B (mT) 

FIG. 7: The surface energy of the liquid layer for increasing 
(open squares) and decreasing (circles) magnetic induction. 
The full circles mark the increase of E s through the successive 
generation of solitons at B = 8.91mT (see also inset). To 
reduce the influence of the perimeter spikes, only the area 
r < 0.88 R was covered, where R is the radius of the vessel. 

Could solitons serve as building blocks in the forma- 
tion of periodic patterns? We have observed that, if addi- 
tional care is taken to suppress the edge-induced inhomo- 
geneity of the magnetization, solitons can form molecule- 
like clusters (Fig.JHl) This may seem to contradict the re- 
pulsive nature of the dipole-dipole interaction; however, 
there is a simple mechanism that can account for the 
binding. Indeed, each soliton is surrounded by concen- 
tric dips representing ring-like regions of depleted mag- 



netic induction. The innermost, deepest, ring is clearly 
visible in Figs. l4l5( and El a higher-resolution measure- 
ment allows to discern another, shallower dip of larger 
radius. The dips create a potential relief which may cap- 
ture the partner soliton(s). As the density of solitons 
grows, the multisoliton cluster evolves towards the hexag- 
onal Rosensweig lattice. It still remains to be understood 
whether the cluster-lattice transition requires additional 
excitation energy. 




FIG. 8: The solitons can form molecule- like clusters: (a) di-, 
(b) tri-, (c) tetra-, (d) penta-, (e) hepta-, and (f) oligomers. 
The height is indicated by switching between black and white 
after each mm. Here B — 8.71mT; each panel covers the area 
of (87mm) 2 . 

In conclusion, we found stable 2D solitons on the 
surface of a ferrofluid in the hysteretic regime of the 
Rosensweig instability. These objects are easy to gener- 
ate and control and they are easily set in motion; this 
opens ways for studying their binding and scattering. 
Due to the conservative nature of the ferrosolitons, and 
unlike the localized structures observed previously in dis- 
sipative systems, the balance of dissipation and energy 
gain plays no role in their stabilization. Instead, we sug- 
gest a stabilization mechanism which appeals to the lock- 
ing of the soliton to the wavelength imposed by the non- 
monotonic dispersion relation. This mechanism can also 
be at work in other conservative systems with preferred 
wavelengths, e.g. in electrostatics and elasticity [l3|. 
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